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Abstract 

Eh 1 

We study the problem of representation of statistical data (of any origin) by a complex 
' probability amplitude. This paper is devoted to representation of data collected from measure- 

Q" 1 ments of two trichotomous observables. The complexity of the problem eventually increased 

comparing to with the case of dichotomous observables. We see that only special statistical 
data (satisfying s number of nonlinear constraints) have the quantum-like representation. 

^ ; 1 Introduction 

The problem of inter-relation between classical and quantum probabilistic data was discussed in 



8, 14|[T5]. We are interested in the 
by complex probability amplitude, 



, numerous papers (from various points of view), see, e.g., [T 

problem of representation of probabilistic data of any origi 
so to say a "wave function". This problem was discussed in very detail in |17j . It has two sources. 
One is purely mathematical: to describe a class of data which permits the quantum-like (QL) 
representation. Another reason to create the QL-representation of data is not so straightforward 
as the first one. In [18] A. Khrennikov presented a hypothesis that biological systems might use 
complex probabilistic amplitudes ("mental wave functions") in processing of statistical data. If 
this hypothesis is correct, then these amplitudes can be reconstructed on the basis of collected 
experimental data. In psychology this approach got the name "constructive wave function ap- 
proach" . 

A general QL-representation algorithm (QLRA) was presented in [T7]. This algorithm is based 
on the formula of total probability with interference term - a disturbance of the standard formula 
of total probability. Starting with experimental probabilistic data, QLRA produces a complex 
probability amplitude such that probability can be reconstructed by using Born's rule. 

Although the formal scheme of QLRA works for multi- valued observables of an arbitrary dimen- 
sion, the description of the class of probabilistic data which can be transfered into QL-amplitudes 
(the domain of application of QLRA) depends very much on the dimension. In |19) the simplest 
case of data generated by dichotomous observables was studied. In this paper we study trichoto- 
mous observables. The complexity of the problem increases incredibly comparing with the two 
dimensional case. 

Finally, we remark that our study is closely related to the triple slit interference experiment 
and Sorkin's equality [16] . This experiment provides an important test of foundations of QM. 



* Research fellowship of Swedish Institute. 

lr Lnus it need not be produced by quantum measurements; it can be collected in e.g. psychology, see |18| . 
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The scheme of presentation is the following one. We start with observables given by QM and 
derive constraints on phases which are necessary and sufficient for the QL-representation. Then 
we use these constraints to produce complex amplitudes from data (of any origin); some examples, 
including numerical, are given. 



2 Trichotomous incompatible quantum observables 
2.1 Probabilities 

Let a and b be two self-adjoint operators in three dimensional complex Hilbert space representing 
two trichotomous incompatible observables a and b. They take values a — ai, i — 1,2,3 and 
b = 0i,l = 1,2,3 - spectra of operators. We assume that the operators have nondegenerate 
spectra, i.e., a, ^ ^ ^ j- Consider corresponding eigenvectors: 

ae a a . =aie a a ., be% = Pie%. 

Denote by = \e a ai ) (e a ai \ and F» t = \e%)( e h l3i | one dimensional projection operators and by P£ 
and Pp the observables repressed by there projections. Consider also projections 

(1) 
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Here the observable P£. = 1 if the result of the a-measurement is a = ai and P£. = if a ^ ai. 
The observables P p are defined in the same way. We have the following relation between events 
corresponding to measurements [P£ = 0] = [P* 2 = 1] V [P^ = 1], [P£ 2 = 0] = [P£ = 1] V [P° 3 = 
!]. [ p a 3 = 0] = [-Pai = !] v [ p a 2 = ■ Thcre are S iven ( b y the QM-formalism) the probabilities 

P b ^PT P (b = f3) = \\P b ^\\ 2 = Me b p )\ 2 , (2) 
p a a ^P^a = a) = \\P«n 2 = Me a a )\ 2 - 
There arc also given (by the QM-formalism) conditional (transition) probabilities 

p% ee P^(b = P\P« = 1) = HP^VIIVII^II 2 = I(e^l4>| 2 . (3) 

We remark that non degeneration of the spectra implies that they do not depend on ip. More- 
over, the matrix of transition probabilities is doubly stochastic. There are also given (V'-depend) 
probabilities 

II p b (p a _|_ p a )ih \ \ 2 

PL 3 - W - MPS* = 1 V f» =l)=P+{b = = <>) = {lpa\paj^ ■ 

where ij, k £ {1, 2, 3}, j, fc 7^ i. We have 

II +^J^H 2 _ II \e b p )(e b p \e%.)(e%.\^) + |e»)(e*|e£ fc )(e£» || 2 

IK^+^JV-II 2 \\\e a aj ){e%^) + \e% k ){e% k \^)\\ 2 

_ ll|e^)|| 2 |(4K J >(eSjV') + (e^e^)«JV')| 2 _ \{e b p \e%.){e%. |V> + (e||e« fc )(e« fc |V>| 2 



KeSJ^)l 2 + l(eSj^)l 2 KeSJV>l 2 + l<eSjV>l s 



Note that p^ 1 
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2.2 Probability amplitudes 

Set ipp = (ip\e b p). Then by Born's rule 

p b P = m 2 - (4) 

We have 

v = ]T#4. ( 5 ) 

Thus these amplitudes give a possibility to reconstruct the state. We remark that -0 = J2 a Pa^i 
and, hence, 

(6) 

a 

Each amplitude V/3 can be represented as the sum of three subamplitudcs 

a 

given by 

^ Q = (^°M> = (il>K)(e a a \e b p ). (8) 
Hence, one can reconstruct the state ip on the basis of nine amplitudes ipp a - We remark that 
l#,aj 2 = K^le^XeSje*,)! 2 =0^. In this notations 

p£„=l^+^ ( a,| 2 /(PS,+PSJ- (9) 



Here = \Jp% i p'^ ai and therefore 

where <~Pji lOLi = asgijjp iai . Moreover, 

(VKJ = v^VK = \[ptZe i6p ^- (10) 

Hence, 

¥?fta m = ^,o m + (11) 

We have a system of equations for phases ip0 iai for i,j,k,l e {1, 2, 3}, 



+ 

We set 



IVa^+^I 2 = \yJpsA^e**°* + sJtfjfU^ I 2 ( 12 ) 

Xljj = COs(ipp iai - ip[3 iaj ) (13) 

and we have 

x (P^+P Q aj Qpt,-(^Pt.+P^Pt 3 ) 
We call A;.ij for the coefficients of interference. 
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2.3 Formula of total probability with interference term 

By using the definition of the amplitude ^/3 iai — (V'l e S)( e al e ^) we obtain 

P% = l^+^+^a*! 3 = l(^|eS s )(eSJe^) + (V|e^)(eS 3 |e^) + ^|e^)(e^|e^)| 2 . (15) 

= mc,)KM\ 2 + mc 3 )KMD\ 2 + mc k )KM\ 2 + 

+ (^l< 3 )(e«j4 ! )(4J<J(CI^ + (<JV')(4 i l< l )(^|e^.)(e^|e^) 
+<<jV'><e^e« i )^|e«J<e»Je^ + <e»J^<e^e«^<V;|e^}<e«Je^). 
Finally, we obtain 

b a b\a , a b\a . a b\a G \ 

Pfl, = PaMcn + PafPptoi + P" h Pf3 iak ( 16 ) 

+ 2cos(^ fta< - Vpi aj )\/p a ,P a] p b fr at P b fr a] +2cos(^ fta!j -%«J\/fewWl 

+ 2cos(p ftaj - <^/3, Qfe ) ^ P aj p ak p b fr a] P b fr ak ■ 

This is nothing else than the formula of total probability with the interference term. It can be 
considered [21] as a perturbation of the classical formula of total probability 

b a b\a . a b\a . a b\a 

= f&Pfta, + PaiPfca, + *W/J,a* ■ ( 17 ) 

If all coefficients of interferes Aj.y = 0, then ([T6|l coincides with (fl7|) . 

2.4 Sorkin's equality in conditional probabilistic form 

We will derive Sorkin's equality by putting (fT3|) in (fl6|) . 

A = p a ai p b L + piA 3 + p a a yt k + 0«« + ^ )#a iv , - (p^ pfe, + )) 

= (pS 4 +*S>£ mvj + ba, + P«>!fe iv)s + «■ +p«>aU 

+ Pa k \Pp,a ivk + Pft QjVi PftaJ- 

This gives us the following constraint on the probabilities 

p%=pUp b L 1 +p b L k -p b L) ( 19 ) 

+ n Q (T) b|a 4- T) b|Q - n b|a *l + n a fr) b|a -I- n b|a - n b|a 1 
+Pa i tf>,9 I a 4 . f + Pf) iajh Pp iaj ) +P ak {Pf) iazk +Pp iajk PftaJ- 

This equation coupling various quantum probabilities can be considered as incrypting of Born's rule 
by using the language of probabilities. This is the discrete version of famous Sorkin equality |22I23] . 
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2.5 Unitarity of transition operator 

We now remark that bases consisting of a- and 6-eigenvectors are orthogonal; hence the operator of 
transition from one basis to another is unitarity. We can always select the 6-basis in the canonical 
way 

/ 1 \ / \ / \ 

(20) 




In this system of coordinates the a-basis has the form 



e a ai = ( Vm^ el6ff2ai > <S a = Vm^ ete " 2a2 > e S 3 = Vm^e ie ^ j (21) 

The matrix 

/ Vm^e i9 ^r ^/p^e* *** ' 

U = y/m^e te ^i y/pp^e 1 ^ Jpp^e**-* 

\ V^W^ftxn y/p^e ie ^2 y/p^;e ie ^2 

is unitary. Hence, we have the system of equations 

E VP/^P/W i(9pmQi ~ <W ° 3 ' ) = ( 22 ) 

or 

2J \fVl3 m a z Pl3 maj COs(6p mCH - 9f3 maj ) = 0, (23) 
m 

^ y/PPnCLiPPmOtj ^(9/3^ ~ 6 ] m<Xj ) = 0, (24) 

m 

where (i 7^ j) (For i = j, the unitarity condition is equivalent to normalization of the sum of 
probabilities by one) +. We now recall that the phases of the basis vectors e° are coupled with 
the phases of the amplitudes ipp a by Hence, we obtain a system of constraints on the later 

phases 

\JVfi mai PH ma] COs[(^ m „, - <j)f3 maj ) + {ta, -CaJ] = 0, (25) 

m 

VP^atP^a, Bia[((j)p mai - <t>p maj ) + (£a, - CaJ] = °- ( 26 ) 

m 

Thus 

COs(£ Qj . - gajy] y/P/3 mai Pt3 m a 3 COs((j)p mQt - (f>{s m0lj ) (27) 
m 

- sin(^ aj - £ajy^ \fPl3 m a z Pp m a 3 sva((j)p mai - </>/3 mCtj ) = 0; 

m 

COs(£ aj - Cajyi y/P/3 m aiP/3 maj Sm(<f>p mCti - <f>p maj ) (28) 
m 

+ sin(£ Qj ~ £ a J ^ \JPp m a.iPp m a.j ™(^„°i ~ <t>P m aj) = °- 
m 

Suppose now the following equalities hold 

2J \JPP m a t Pf3 maj COs(0/J mQl - 0/3 m a 3 ) = 0, (29) 

m 

^ <s/PP m aiP0 m <Xi Sin(^ m ai - <l>i3m<Xj) = 0- (30) 
m 

Then by {27 } , {25 } the equalities {23]), (J2H) and, hence, {22} hold. Thus conditions {25} , {30} 
imply unitarity of [7. It is clear that in turn {23} . {24} imply {29} . {30} for arbitrary f Q . Thus the 
later conditions are equivalent to unitarity of U. 



5 



3 Mutually unbiased bases 



In previous considerations we introduced the coefficients of interference Xnj on the basis of the 
phases <fp iai , "^ftay by A;.y = cos{ipp iai — ipp iaj ). However, we see that they also could be defined 
by using only probabilistic data, see ([T4"]) . Can we the go other way around and to find phases 
fPiat, PPiaj on the basis of the interference coefficients given by (fT4|) ? We will study this general 
problem in section [5] Now we would like consider an example. To be sure that the problem has 
a solution, we start with data and the corresponding interference coefficients generated by QM. 
In section [S] we shall operate with statistical data an arbitrary origin. First we show that there 
exist probabilistic data such that (|2T))) and (|3U| hold, therefore consider the following situation. 
Let p/3 iaj =1/3 where i, j = 1, 2, 3. We will also put 

0p 2 ai — Offsets — —9/3 2 a2 = —Opsoti = 27r/3, (31) 

where all the other 9/3 iaj = 0, i,j = 1,2,3. We insert this in the basis in (|2"Tj) and obtain the 
orthonormal basis 




where w = e l27T / 3 . Bases (|2U)) and (|3"2")l are mutually unbiased, i.e. ,|(e|.|e£ 



Note that \ip\ 2 = 1. Then by equation @ and after some calculations this gives that, 

Pai = l( e /3j^)| 2 = g (3 - cos (71) - cos (71 - 72) + 2 cos (72) + VSsin^i) + V3sin(7i - 72)) , 

1 _ _ (33) 

Pa 2 = l( e ^ a |V')| 2 = g (3 ~ cos (71) + 2 cos (71 - 72) - cos (72) - sin (71) - V% sin (72)^ , (34) 

Pa 3 = I <e^ 3 1-0> I 2 = g (3 + 2 cos (71) - cos (71 - 72) - cos (72) - \/3sin(7i - 72) + \/3sin(7 2 )) . 



(35) 



G 



It is straightforward to see that p ai + p a2 + p a3 
l,k,j = 1,2,3, k / j are calculated by ([§]) 



P/3iai2 



1. The conditional probabilities pp l0lkj 

4 cos (71) + cos (71 - 72) - 2 cos (72) + a/3 sin (71 - 72) + 2 a/3 sin (72) - 6 
3 (-2 cos (71) + cos (71 - 72) + cos (72) + \/3 sin (71 - 72) - \/3 sin (72) + 6) 
2 (cos (71) + cos (71 - 72) + cos (72) + a/3 sin (71) + a/3 sin (72) + 3) 
3 (cos (71) - 2 cos (71 - 72) + cos (72) + V3 sin (71 ) + V% sin (72) + 6) 



for 



(36) 



Pf3ia 2 3 



■ — 



G 



P/32«13 



'1 - 72) + 4 cos (72) + 2\/3sin (71) - a/3 sin (71 - 72) - 6 
~ ' "(72) - V3sin(7i) 
72) + 2 a/3 sin (71 
(72) + \/3sin (71 - 

j (72) - 2\/3sin (71J -t- y 6 sin ^72 ) - < 

- 2 cos (71 - 72) + cos (72) + a/3 sin (71) + a/3 sin (72) + 6) 



— 2 COS (71) + COS (71 — 72 J T t ™ ^72; T sm \'J1) — \ o sin ^71 — r y2) — u 

sin (71) - \/3sin (71 - 72) + 6) 
4 cos (71) + 2 cos (7! - 72) - cos (72) + 2 a/3 sin (71 - 72) + a/3 sin (72) + ' 
(71) + cos (71 - 72) + cos (72) + \/3 sin (71 - 72) - a/3 sin (72) + ' 
71) + 4 cos (7x - 72) + cos (72) - 2\/3sin (7!) + a/3 sin (72) - 6 

3) 



3 (cos (7; 
2 (cos (71) + cos (71 - 72) + 



- 1. w.> : . [ , t ^ . > 1 - > 1 1- cos (72) - a/3 sin (71) - a/3 sin (71 - 72) t ■>< 
3 (cos (71) + cos (71 - 72) - 2 cos (72) - \/3 sin (71) - a/3 sin (71 - 72) + 6) 
2 (cos (71) + cos (71 - 72) + cos (72) + a/3 sin (71 - 72) - a/3 sin (72) + 3) 

"sin' 



2 (cos (71) + cos (71 - 72) + cos (72) + a/3sL, W1 v , 
3 (-2 cos (71) + cos (71 - 72) + cos (72) + a/3 sin (71 - 72) - a/3 sin (72) + 6) 

_ cos (71) + 4 cos (71 - 72) + a/3 sin (71) - 2 (cos (72) + V3 sin (72) + 3) 

3 (cos (71) - 2 cos (71 - 72) + cos (72) + a/3 sin (71) + a/3 sin (72) + 6) 
cos (71) - 2 cos (71 - 72) + 4 cos (72) - a/3 sin (71) + 2 a/3 sin (71 - 72) - 6 
3 (cos (71) + cos (71 - 72) - 2 cos (72) - a/3 sin (71) - a/3 sin (71 - 72) + 6) 

probabilities are found. Moreover we find that p\ =1/3, I — 1,2,3 by inserting (fT3")) ir 
lirectly by (fT9|) . We let 71 — 72 in order to get more compact expressions. This leads to 



All 

or 



Bia 12 



82012 



83012 



.11 piuutiuiniics aic louiiu. iv 

r directly by (fl9|). We let 71 

4sin (72 + f) -5 
3 (2 sin (72 

-5 cos (72) + a/3 sin (72) + 8 

3(2sin( 72 + f)-7) ' Pf32ai3 = 

2 (-2 cos (72) + \/3 sin (72) - 4) 

3(2sin( 7 2 + f)-7) ' mai3 

Here we calculate \ } jk when 72 = 71 



4 sin (72 



I = 1, 2, 3 by inserting ([T5|) in (fTol) 

(37) 



l)-5 



3(2sin (72 + f^ -7 



1 

6' 



02^23 



1 

6' 



83023 



2 (-2 cos (72) + \/3 sin (72) - 4 
3 (2 sin (72 + |) -7) 
-5 cos (72) + a/3 sin (72) + ■ 



3 (2 sin (72 + I) -7) 



Al,l2 — Al,23 — 



1 + sin ( 72 + f ) 
10 -8 sin (72 + f) 



-^2,12 — A323 — 



4 cos (72) + 2A/3sin (72) + 1 



2^4 sin (f - 2 72 ) + 2 sin ( 7z + f ) + 6 



Al,13 — 1 



A243 — A342 — —7: 



A2,23 — A342 — 



-5 COS 



2, 1 A sin (f - 2 72 ) 



(72) + a/3 sin (72) - 1 



-2sin (72 + f) +6 



2 Recall that p b J a = r> b } a 
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We calculate numerically the extreme values for all A to prove the existence of angles between 
and 2n for arccosA. Thus, calculate ^-A = and the limits when 72 — > ±47r/3, 

-T-Ai,ia = -^-Ai. 2 3 = => 72 = J + 2ttC*, C G Z 
CZ72 (172 ' 3 

and A142 = Ai,23 = —1,72 — f + 27rC. The problems arise when 72 — > Att/3 for A142 and Ai,23, 
since lim 72 ^ ±47r /3 = 0. Let A142 — Ai,23 = when 72 = An/3. Moreover there exist angles 

7!"/2 < <pp uai - Vp u a 3 < 7T) ( 38 ) 

where arccosAi i i2 = fp lt ai ~ <P/3i,a 2 - We see that the denominator of A2,23, Aa^, A3 ; i2, \2,i2 
goes to zero when 72 — > 47r/3±. We therefore examine the limits 

lim A 3 . 23 = lim A2.12 = hm A2.23 = lim A 3 ,i 2 = ±— . (39) 

72— *4ir/3± 72— >47r/3± 72— >4ir/3=F 72— >47r/3=F 2 

Analysis of extreme values provides us with the following 

-^-A 2 ,23 = -^-A 3 ,i2 = =>• 72 - 2ttC7, CeZ, (40) 
d72 072 

and 

^-A 3 . 23 - ^-A 2 ,i2 - 72 = %■ + 2ttC, CeZ. (41) 
«72 "72 3 

The minimum values are given by (|3T))) . minAs^a = minA3_23 = rninA3,23 = minA3,23 = — 72 — > 
47r/3 or 72 — > — 47r/3. The maximum values are given by (|40p and (|4ip . maxA2.23 = maxA3 j i2 = 
1, 72 = 2ttC, C G Z. maxAa^ = maxA2.i2 = 1, 72 = + 2ttC;C € Z. This prove that there 
exist angles such that 

(42) 



< 




- m 


.03 


< 


5tt 

IT 


< 




- v% 


,C*3 


< 


IT 


< 




- ¥>/3 2 


,C* 2 


< 


IT 


< 






. a 2 


< 


IT 



and 



J 1 ,o,=0, (43) 
2tt 



2tt 

^3,01 - <P03,a2 - -y 



4 Triple slit experiment 

An interesting example of interplay of two incompatible trichotomous observables is given by the 
triple slit experiment - a natural generalization of the well know two slit experiment. There are 
given: a) a source of quantum systems which has very low intensity (so it might be interpreted 
as single-particle source); 2) a screen with three slits a = 1, 2, 3 and each of them can be open or 
close on the demand; c) registration screen; typically it is covered by photo-emulsion; this produces 
the continuous interefrence picture; we shall consider discrete experiment. The a-observable gives 
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slit's which is so to say is passed by a particle on the way from the source to the registration 
screen. To measure a, an experimenter puts three detectors directly behind slits. We set a = oq, 
if the detector behind the ith slit produces a click. By the assumption the source has so low 
intensity that one can neglect by double clicks (e.g., the detectors never click simultaneously). We 
can find probabilities p^ot = f,2,3. This is the first experiments producing a-probabilities. Now 
we consider basic experiments. 

To make the second variable discrete, we put detectors in three fixed places of the registration 
slit. It gives us the observable b = 0,0 = 1, 2, 3. Thus = 1 if the first detector clicks and so on. 
The experiment is repeated at a few incompatible contexts, see [17] for general presentation 

C123 : all slits are open; probabilities of b detection collected in this context are probabilities 
P b p = Ptl>{b = P) from section I2TT1 In the QM-formalism context C123 is represented by a quantum 
state ip. 

C a ,a — 1,2,3 : only the slit a is open; corresponding probabilities axe p b g^, transition proba- 
bilities. In the QM-formalism context C a is represented by the quantum state e°. 

C aij ,i 7^ j ■ only slits a, and otj are open (so the slit a>k,k ^ is closed); corresponding 

probabilities are PajP^a ■ I R the QM-formalism context C aij , i 7^ j is represented by the quantum 

state ^ = n£;t£;£ir 

Thus all probabilities discussed in section |2~T1 can be collected in this experiment. It is possible 
to check whether these experimental probabilities match the predictions of QM. The easiest way 
for experimenters is to check Sorkin equality (1191) . 

Recently the group of Gregor Weihs performed the triple slit experiment, see [16] . They cliam 
that Sorkin's equality and Born's rule are violated by their experimental statistical data. 

5 Construction of a complex probability amplitude satisfy- 
ing Born's rule 

Now we have a pair of trichotomous observables a and b taking values a = a,-, i = 1,2,3 and 

6 = 0i, I = 1, 2, 3. We do not assume that they have any relation to quantum physics; e.g., these 
are some random variables observed in biology or finances. It is assumed that there are given 
probability distributions of these variables 

p* i =P(b = l ), p a ai =P(a = ai ). 

Thus 

3 3 

E^ = 1 'E^ = L ( 44 ) 



1=1 1=1 

It is also assumed that there are given conditional probabilities Pg* a . = P(b = fii\a = a 2 ). We 
know that for any sort of data the matrix of transition probabilities is stochastic, i.e., for each 

I>£n=1- (45) 
1=1 

Finally, we assume a possibility to collect the data on measurements of observables P£. , i = 1, 2, 3, 

probabilities Pp" ak . = P{b — 0l\Pa t — 0)- "The detector corresponding to the value a = on does 
not click, so the value of a is either a = aj or a = ak, where j, k ^ i. However, we do not know 
the value of a. In this context we measure the b- variable." Fo any sort of data, we have 

3 



E*« =1 - (46) 



1=1 



3 It is surprising that it has not been done for long ago! 



9 



5.1 Complex amplitude matching Born's rule for one observable 

Now we want to find a complex probability amplitude ipp t such that Born's rule (for the b- variable) 
holds: l^g, | 2 — P b p r We copy the QM-scheme, so we represent ipp, — ip/3 iai +i>p ia2 +*Ppia3 j where the 

sub-amplitudes V'ftai = \/P'i i P b faa- elipi3iai an d phases are determined by the system of equations 



(I13P . It is convenient to work with the interference coefficients, see |17j , given by right-hand sides 
of these equations 

, = _ (^+K J >t,-(^t.+K J pt 3 ) U7) 

Interference coefficients obtained in quantum physics are always bounded by 1: 

\\ij\<l. (48) 

However, since we start with data of any origin, the condition (|48fl has to be checked to proceed 
to representation of data by complex amplitudes^] If the system of equations, m = 1, 2, 3, 

(49) 









1*^2) 


— A mj i2 










= Kn : 23 


COs(<^ /3 „ 






lOs) 


= A r „43. 



has a solution (three phases) then we can construct the probability amplitudes ip/3 iai and, hence, 
the probability amplitudes tj}p l and the corresponding vector ip. 

However, in general such amplitudes will not provide a solution of the "inverse Born problem" , 
namely, Born's rule can be violated. To obtain the real solution one should solve the system (|49l) 
under the constriant (1191) . Thus, to proceed toward a proper complex amplitude, one should first 
check the validity of (fl9|) and then to solve the system (|49[) . It is convenient to express "triple 
probabilities" P/3 m , aij through coefficients of interference 

PPmaij = T {Pa % P0 m a, + Pa.jPp m a j ~ 2 A^ may s /Va. i Pp m a i Pa. j P/3 maj ) • (50) 

Pcti ~r Polj 

We remark that if (|4"8")l holds, then triple probabilities given by (IBTfl) are always nonnegative. By 
using the A- variables normalization equations (|46[) can be written as (j, k = 1, 2, 3) 

3 

^ *l,jky/PP, aj Pf} tak = 0. (51) 

z=i 

We also can write Sorkin's equality (in fact, the formula of total probability with interference 
terms) as 

p% = piA, + p^Ai + p a a A k + 2A '.« \Ip^AA 3 ( 52 ) 

+2A Wfe y PaiP ah p^ ai P^ ah + 2 *hik\jp a] p ak p 0iao P 0iak - 

Hence, to obtain Born's rule for the ^-variable which matches the intereference formula of total 
probability, we have to find A satisfying equations (|51j) and (|52j) and put such A into equations 
(|49| . then solve this system of equations. In general, it is a complex problem. 
Thus, finally, we can write the complete system of equations: 

3 

x Ljk^/Pp ia] Pp iak = 0, j, k = 1, 2, 3; (53) 



4 If this condition is violated then data may be represented by so called hyperbolic probabilistic amplitudes I2UI . 
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iC«l 


- <pp„ 


1*^2 ) 


= ^m,12 


cos (¥>/9„ 




- <p/9„ 


i«3 ) 


= A Ifl! 23 


cos^ 


I Oil 


- Pft, 


.03 ) 


= A m! i3 



= p a a /L + PlA 3 + PlA k + 2A w yJpoiPoAA,; (54) 

i o \ / ''I' 1 ''I' 1 i <"> \ / b\a b\a 

(55) 
(56) 
(57) 

Solution of this system will provide us a complex probability amplitude ip such that | (ip\eg) \ 2 = 

!>'•■■ 

Let us consider the case of maximally unbiased matrix of transition probabilities; 

pf w =V3, V/,j = 1,2,3 (58) 
Moreover, to simplify the task by the factor of three, we will put all 

p% = 1/3, VZ = 1,2,3 (59) 
Now, let us introduce new variables x > and y > 0: 

\Jp%Jp% 2 = x > \Jp%Jp% s = Vi \Jp%Jp% 3 = y/ x , ( 60 ) 

That means that 

2 2 2 2 

_a _ X V „a V a ^ , fi1 x 

"a 2 „,2™2 i „,2 i „.2 > „,2„2 i „,2 i „,2 ' \ 01 / 



Q ' y 2 x 2 + a; 2 + y 2 ' Q2 y 2 x 2 + x 2 + y 2 ' Q3 ?/ 2 a; 2 + a;^ + y l 

and the condition + + = 1 always holds. Let proceed for a particular choice of 
interference coefficients (ansatz) A/^2 = [i, \i,v& = — /!,, thus A/, 23 = 1 — 2/i 2 by 



A ; , 23 = A U 2A /il3 ± ^(1 - A 2 12 )(l - A 2 13 ). (62) 
The system of equations (|47)) for A under conditions (|58t and (l59t have the form: 



A,i3 = f (y+i) (pfc I3 -^) 

\ 3 /y x\ / b | a 1 
A ^ 23 "2 U l>*-~3 



We write this as: 



b\a _1 2^ f,| 1_ 2/i b | a 1 2 (l - 2^ 

9 Aa 12 3 3 (rc + i) ' "3 3 ( y + 1 V ~ 3 + » ^ ■ «\ ' 



3 (f- 









where /i is a parameter. The probabilities given by satisfy the relation in (IT91 which in this 
case looks as: 

2 

(a; 2 + l)y 2 pp iai2 + (y 2 + l)x 2 p [jiai3 + (x 2 + y 2 )p/3 ia23 = -{x y + V 2 + x2 ) (65) 
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{6.925, 2.225} 




Figure 1: Left side of the expression (|68|) is plotted versus parameter v = A;, 13 for the fixed values 
of x = 6.925, y — 2.225 and different signs before square roots in expressions for A/.23, see (|6T)|) , and 
A;. 12, see (|70|) . We can see, that we can often satisfy (|68|) by properly choosing signs in expressions 

m and rm. 



Putting into (|rJ5"j) . we get an equation for fi: 



o 2 , / \ , n _^ (y - a;) ± V - y) 2 + 8 ,.„, 

2/i + (x — — 1 = => (i = (66) 

We are interested in the case then all absolute value of lambdas are less than one 

|A Mj | <l,Vy = 12, 13,23. (67) 

It is satisfied when \fi\ < 1. So, for the case of \x — y\ < 1 , for both roots of (|r?6")) conditions 
I Ai,ij I < 1 are valid, if x — y > 1, then (|66p with the plus sign suits \Xi t ij \ < 1, otherwise y — x > 1, 
then (|66[) with the minus sign is valid. 

Now we proceed in a general case, i.e. without ansatz A;. 12 = /i, A;. 13 = — /1, A;. 23 = 1 — 2/i 2 . 
Conditions ([58| and (|59|) equation (fl9|) . which is equivalent to Born's rule, comes down to: 



yXl,12 + X\l t l3 + A;,23 = (68) 

We should combine it with the constraint, see (|6"2")l for A;. 12, A;,^, A;. 23 to have simultaneous solu- 
tion 



Ai, 23 = Ai.12Ai.13 ± - Af jl2 )(l - A 2 13 ) (69) 

We have two equations for three variables, thus we can express the solution as a one-parametric 
family. Let us choose v = Aj.i3 as a parameter. Then 



-xv(y + v)± yj (v 2 - l)(v 2 x 2 - y 2 - 1 - 2yv) 

A U2 = 2 , (7U) 

y + 1 + 2yw 

and A;.23 can be obtained from equation (|62j) . We have to make sure that A;. 12, A;,i3, A;,23 exist, 
are real and satisfy (|67|) . given real and positive x and y. In this, more general, case 

b\a _ 1 , 2A;.i2 f,| a _1 2A;.i 3 b \ a _1 2A;, 2 3 /„-. v 

~ 3 + 3 (x + i) ' PftQl " 3 + o („ , 1 V PftQl ~ 3 + q f * m 1L\ ' 1 j 



Seeing that all values in the parentheses in (|64p are greater than 2, each of the this probabilities 
non-negative and smaller than 2/3, if Az.12, A/.13, A; ; 23 are smaller than 1. The main problem is 
to describe possible rangers of parameters in (|70j) which give us < 1, see figure 1-3. W e 

remark that A;, 13 is a a parameter, \v\ < 1. 
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{0.555} 




Figure 2: Left side of the expression (|68[) is plotted versus parameters cc, y, for the fixed v = Xi,ia = 
0,555. Minus sign is taken before square root in (f70|) . plus sign is taken in (|69l) when obtaining 
Az,i3 values. We can see that it is equal to zero as (|68j) demands on a large scope of x, y parameters 
values. 



5.2 Complex amplitude matching Born's rule for two observables 

We now want to select an orthonormal basis e a such that, for the state ip constructed in the 
previous section, K^le^)! 2 — P a - We turn to considerations of section l2~5l Since vectors of this 
basis can be selected up to e^=. We can select 6p iak — ( f>i3 1 a k , i.e., set £ Q = 0. Of course, to 
guarantee orthogonality of this basis, constraints (l29l) . (|30f should be taken into account: 

^ ^PPmaiPPrr.a, X m ,ij = 0, (72) 



e m,ijy/Pfi mai Pp maj ^J 1 - A^y = 0, (73) 

m 

where e m ,ij are signs which are selected in a proper way. 

Moreover, the matrix of transition probabilities has to be doubly stochastic, i.e., besides (|45|) . 
we should have 

X>&, = 1 ( 74 ) 

1=3 

for each i = 1,2, 3. 

Under these conditions the complex amplitude ip produced by our algorithm matches with 
Born's rule for both obsrevables, a and b. 

Example 1. We can take 

Al,12 = ^1,23 = — M; Al,13 = 0; 

A2,12 = —fit ^2,23 = 0, A2,13 = fJi] 
A3, 12 = 0, A3, 23 = Pi A343 = — fi. 
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Figure 3: Value of \i,i2 as function (f70|) of parameters x,y, for the fixed v = Xi,i3 = 0.585. Plus 
sign is taken before square root in [70]). One may see that value A;,i2 is within (—1, 1) limits for 
the most part of the values x and y. 

This gives us the solution /i = First take the plus-case(e.i. fj, = i^j)- We select tfi m x = 
v mi m = 1,2,3. For /3%, we obtain the system of equations: cos(z/i — <j>i2) = 7^jCOs(^>i2 — ^13) = 
_ 7f' cos (^i _< ^i3) = 0- Hence, <j>u = v\, and </>i 2 = vi— 7r/4, </>i 3 = z/i+7r/2 or i2 = ^i+7r/4, 0i 3 = 

!/ X - 7T/2. 

Then, for^ 2 , cos(^ 2 -^22) = -7751 cos(<?!> 22 -</> 23 ) = 0, cos(^ 2 -</> 23 ) = Hence, 4> 2 i = V2, and 
</>22 = z>2 + 37r/4, 023 = ^ 2 +7r/4 or 22 = f 2 - 37r/4, 23 = 2/ 2 -7r/4. Finally, for /3 3 , cos(t' 3 -^ 32 ) = 
0, cos((/) 32 -0 33 ) = ^,cos(iAj-0 33 ) = -^=. Hence, </> 3 i = i/ 3 , and 32 = i^ 3 + 7r/2,0 33 = v 3 +3n/4 
or (j)32 = — it/2, 33 = v-i — 3tt/A. We have from equation JTJ),© and ([TUf that 

where £ Q . = and ■\/p% rn P(3 l a rn = 3- Thus 

^ - ^ + e Wl ^ /4 + e-^'/ 2 ) = ^(1 + e^/ 4 + e 4 ^/ 2 ) = ^((1 + ^) ± i(l - ^)); 

tip, = ^T-i 1 ± *V2) and > nnall y> i>Ps = ^T-d 1 ~ 75) ± *(* + 75))- We remark that |^ ft .| 2 = 
1/3, j = 1,2,3. 

QLRA produces following possible realizations of the "wave function" : 

V = i[((l + -J=) ± i(l - -^))e& + (1 ± zV2)e^ 2 + ((1 - -L) ± i(l + ^))e&], (75) 
where are the orthonormal basis 

e\ = ( Q j , 4 2 = f j , e> 3 = ( ^ j . (76) 

For the minus-case (e.i. fi — — A- ) QLRA produces following "wave function": 

= i[((l - -L) ± i(l + -L))e^ + (1 ± *V2)e£ a + ((1 + -L) ± i(l - ^=))e£,]. (77) 
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